Abstract. Accretion flows onto supergalactic-scale struc- 
tures are accompanied with large spatial scale shock 
waves. These shocks were postulated as possible sources 
of ultra-high energy cosmic rays. The highest particle en- 
ergies were expected for perpendicular shock configura- 
tion in the so-called " Jokipii diffusion limit" , involving 
weakly turbulent conditions in the large-scale magnetic 
field imbedded in the accreting plasma. For such configu- 
ration we discuss the process limiting the highest energy 
that particles can obtain in the first-order Fermi acceler- 
I ation process due to finite shock extensions to the sides, 
along and across the mean magnetic field. Cosmic ray out- 
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' flow along the shock structure can substantially lower the 

■ upper energy limit for conditions considered for super- 
d ' galactic shocks. 
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1. Introduction 

Formation of large-scale extragalactic structures (galaxy 
clusters and super-clusters) involving accretion of plasma 
to the gravitational potential wells is accompanied by for- 
mation of strong shock waves (Ryu et al. 1993). These 
shocks were considered as viable sources of ultra high 
energy cosmic rays (UHECRs), up to energies E„iax ^ 
1018-19 eV (Kang et al. 1997, HiUas 1984). Evaluation 
of this maximum energy is not reliable due to inadequate 
knowledge of physical conditions at the shocks considered. 
The first estimates show that maximal particle energies 
accelerated in such conditions must be limited because 
of the extensive time scales involved. The magnetic fields 
in the considered fiows are not expected to be stronger 
than ^ 1 /iG for the accretion in rich galaxy clusters 
and ^0.1 fiG for supergalaxy scale accretion (Kronberg 
2001). At the same time the former flows can involve 
plasma velocities ~ 1000 km/s, while the latter a few times 
lower ones (Juszkiewicz et al. 2000). The ordinary diffu- 
sive Fermi acceleration considered (Drury 1983) is fastest 
for highly turbulent conditions with the diffusion coeffi- 
cient comparable to the Bohm diffusion coefficient, kb. In 
the considered situations, the involved acceleration time 
scale T ~ '^s/Kiocfc becomes longer than the age of the 
universe for particle energies greater than ~ 10^^ eV. 

However, as pointed out by Kang et al. (1997; cf. 
Jokipii 1987 and Ostrowski 1988), even faster acceleration 
can take place at perpendicular shocks occurring when 
large-scale magnetic field structures are present in the in- 
tergalactic medium and the accretion proceeds perpendic- 
ular to the field direction. Then, the cosmic ray cross-field 
diffusion coefficient governs the first order Fermi accelera- 
tion process at the shock. With the involved spatial scales 
along the shock normal of respectively 1 Mpc for a clus- 
ter, or a few Mpc for a super-cluster flow, there is enough 
space and time to accelerate particles up to ^ 10^^ eV, 
three orders of magnitude greater than in the Bohm dif- 
fusion limit (cf. Siemieniec-Ozi§blo & Ostrowski 2000). 
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Fig. 1. Qualitative sketch of the physical situation con- 
sidered near the extragalactic-scale accretion shock: lines 
represent the mean magnetic field which is (quasi) perpen- 
dicular to the plasma flow at the shock. A non-uniformity 
of the flow results in a finite shock extension along the 
mean magnetic field lS.Xmax- Particle escape along the 
magnetic field from the acceleration region at the shock 
is represented by the introduced free escape boundaries 
('feb') indicated by dotted lines). 

Below, in section 2, we present evaluations of the accel- 
eration time scales in perpendicular shocks. The derived 
time scale is used in section 3 to discuss the particle diffu- 
sion and drifts along the shock front, perpendicular to the 
plasma flow direction. We show that particle escape to the 
sides of the shock structure may substantially limit E^ax 
derived for perpendicular shocks. In section 4 we briefly 
discuss the importance of this constraint for large-scale 
extragalactic shocks. 
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2. Particle acceleration at perpendicular shock 

waves 

The acceleration time scale derived within the diffusion 
approximation for oblique shocks (Ostrowski 1988) be- 
comes formally invalid for perpendicular shocks which are 
always superluminal. Nevertheless, Jokipii (1987) used the 
diffusion equation to derive the energetic particle acceler- 
ation time scale for such shocks: 



Ui-U2\Ui u, 



K2 



(2.1) 



where ki and Ui {k2 and U2) are the normal diffusion co- 
efficient and shock velocity as measured upstream (down- 
stream) of the shock. With the scaling assumed by Jokipii 
(1987) 



Kj_ 



1 



1 + 



(2.2) 



where symbols _L and 1| indicate quantities perpendicu- 
lar and parallel to the mean magnetic field, respectively, 
and A|| is a mean free path for particle scattering along 
the magnetic field (k|| — ^X^^w, where w is a particle ve- 
locity), and rg = pc/eB is the particle gyroradius. The 
equation (2.2) describes the cross field diffusion, where, 
at the mean free time, the particle moves on average A[| 
along the magnetic field and is scattered at a distance Vg 
across it. Continuing the acceleration process at perpen- 
dicular shocks requires sufficiently effective particle cross- 
field diffusion. With r] = X\\/rg this requirement can be 
expressed (Jokipii 1987) as 



w 



3. Maximum energy of accelerated particles 



(2.3) 



3.1. t: a diffusive approach 

In order to evaluate the maximum energy attainable in 
the acceleration process acting at the perpendicular shock 
wave of a finite extension AXmax along the magnetic field 
one has to consider particle escape from the shock region, 
represented by free escape boundaries of figure 1. Let us 
assume Eq. (2.1) to be a valid acceleration time estimate. 
Additionally, we consider conditions with the constant r] 
upstream and downstream of the shock i.e. t] = Ay^i/rg^i = 
A||,2/''g,2 • With the shock compression ratio r = U1/U2 
and a weakly perturbed upstream magnetic field Bi, the 
downstream magnetic field B2 = rBi and the particle 
gyroradius becomes rg,2 = i'g,i/r . Fast acceleration at 
perpendicular shocks appears only for weakly turbulent 
conditions with Ay ^ Tg. Prom Eq. (2.2), in this limit 

K_L = Ky . For simplicity let us consider a purely 

perpendicular shock with the mean magnetic field par- 
allel to the shock surface and, respectively, Ki = and 



K,2 = K_L_2j ^iiid rclativistic cosmic ray particles with w = c. 
Then, with the Bohm diffusion coefficient kb = jJ'gC, the 
diffusion coefficients 



Ki = -kba {i = 1, 2) 
and the equation (2.1) transforms to 



T = 



2r lr„,i 



(3.1) 



(3.2) 



r — 1 f] c \Ui 

With the constraint (2.3) and the strong shock with r = 4 
the acceleration time scale must be longer than c/Ui times 
a particle gyroperiod. 

3.2. T : a lower limit for the perpendicular shock 

An independent lower limit for the acceleration time scale 
can be derived by considering physical conditions ap- 
proaching a weak scattering limit. Then, in the absence 
of scattering, let us consider the shock-drift acceleration 
process of upstream particles at the perpendicular shock 
propagating in a uniform magnetic field. A particle with 
the perpendicular momentum component p± compressed 
in the shock will increase it to 

P±,2=P±,iVr (3.3) 
due to the approximate magnetic momentum conserva- 
tion. For strong shocks, on average particles roughly dou- 
ble their momenta at individual transmissions through the 
shock, which take approximately 



At = 



2r„,i 2r 



5,2 



(3.4) 



Ui U2 

Thus the characteristic acceleration time cannot be 
shorter than r ~ At. The same result can be obtained 
by considering particle drift along the shock, perpendic- 
ular to the weakly perturbed magnetic field. Addition of 
the diffusive process that is required to allow downstream 
particles to return to the shock to continue acceleration 
can only make t longer. One may note that the minimum 
acceleration time in the Jokipii (1987) derivation approxi- 
mately coincides with the above one. Thus, the time (3.4) 
can be considered to be an absolute lower limit for the 
acceleration time in the perpendicular shock and thus it 
defines the ultimate limit for k± reduction (or the upper 
limit for the considered 77) in derivation of the expression 
(3.2). 

3.3. New constraints for E^ax 

Cosmic ray removal from the acceleration region near the 

shock proceeds due to particle advection with the plasma 
ffow downstream of the shock and due to diffusion or drifts 
towards the escape boundaries around the shock region. 
In the present paper we discuss the role of such escape in 
limiting the maximum energy of accelerated particles. 
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3.3.1. A limit due to escape along the magnetic field 

First let us discuss the diffusive escape along the magnetic 
field (cf. Fig. 1) limiting the maximum particle energy at 
the shock of a finite extension. To evaluate this limit we 
consider a solution of the 1-D spatial diffusion equation 
for particles moving along the magnetic field (along the 
a;-axis of our coordinate frame), with the diffusion coeffi- 
cient For this evaluation we neglect particle transport 
perpendicular to this field, both along the shock normal 
and to the sides. For No particles injected at a; = and 
t = Q, their distribution along the a;-axis is given by 



With the above expressions for the perpendicular shock 
the equation (3.8) can be rewritten as 



U L 



c ra 



(3.9) 



In order to double the original energy in this process, 

Ai? = i?, an average particle has to move a distance 
L = (c/U) ■ Tg along the shock surface. Thus, if the shock 
extension available for drifting particles L = AYmax/^, 
the maximum energy of particles accelerated by the shock 
drift mechanism is provided by the condition 



n{x, t) 



^47rK|| t 



exp 



4k;ii t 



(3.5) 



larger 



r9< 



(3.10) 



analogous to the above condition (3.7). 



If a dispersion of this distribution y^2K|| t beco 
than the shock sideway extension AX^axi the accelerated 4_ Conclusions 
particles tend to efficiently escape to the sides, leading 
to the spectrum steepening and, finally, to formation of 
the energy cut-off. Such modification of the particle spec- 
trum is possible when the time required for escape be- 
comes comparable to the acceleration time scale r, i.e. 
when the particle dispersion along the x-axis over time r 
becomes comparable to the shock extension AX^ax' 



(3.6) 



Applying the above Eq. (3.2) and the expression for K|| , we 
obtain the requirement for diffusively accelerated particles 
at such finite-extent shocks as 



AXr, 



U 



(3.7) 



As the particle gyroradius Vg = pc/eB, this expression 
provides the upper energy limit for accelerated particles 
due to particle escape along the mean magnetic field. For 
nonrelativistic shocks with the factor U/c <^ 1, Eq. (3.7) 
provides a severe constraint for the maximum accelerated 
particle energy. 

3.3.2. A limit due to drift across the magnetic field 

As discussed in section 3.2, particle acceleration at the 
perpendicular shock with the weakly perturbed magnetic 
field involves particle drift along the shock surface, across 
the magnetic field. Thus, the shock extension along this 
direction (along the y-axis of our coordinate frame) also 
provides a limit for the maximum energy of accelerated 
particles. The energy gain of a particle drifting along the 
shock is 



AE = \q£L\ 



(3.6 



where the "VxB" electric field £ = -{Ui x Bi)/c = 
— (t/2 X B2)/c, q is the particle electric charge and L is 
the distance the particle drifted along the shock {AE > 0). 



The maximum energy a particle can gain in the drift - dif- 
fusion acceleration expressed by the famous Hillas (1984) 
formula needs more specification concerning the charac- 
teristic spatial scales of the shock. Apart from the shock 
normal extension (Bhattacharjee & Sigl 2000) one also 
needs information about its sideway extent {AX^ax and 
AYmax) influencing the maximum energy and, at highest 
energies, the spectral shape of the accelerated particles. 
As an illustration let us consider an accretion flow with 
velocity U = 300 km/s, the size of /i w 3 Mpc along the 
flow and of L ~ 10 Mpc to the sides, harbouring a shock 
wave inside with, say, AXmax and/or AYmax ~ i . In the 
perpendicular shock configuration the plasma normal ex- 
tent provides the upper limit for particle energy expressed 
as rg < h. On the other hand particle escape to the sides 
provide the limits (3.7) and (3.10). If the magnetic field in 
the plasma equals 0.1 /xG, the former constraint requires 
approximately that Emax < 10^*^ eV, while the others pro- 
vide the limit Emax ^ 10"'^^ eV. When varying the above 
accretion flow parameters within the reasonable ranges 
one finds that the upper energy limit for particles accel- 
erated at such shocks can hardly exceed a scale of EeV, 
and in most cases is close to 10^^ eV. This low E^ax re- 
sults from the fact that for the weakly perturbed magnetic 
field required for fast acceleration at the considered per- 
pendicular shock, the accelerated particles can diffuse and 
drift far along the finite shock to escape to its sides. For 
highly turbulent conditions and/or the parallel shock con- 
figuration a similar value for E^ax arises due to the slow 
acceleration constrained by radiation losses (cf. Kang et 
al. 1997). 

One should note that obtaining the same expressions 
for constraints (3.7) and (3.10) is not a coincidence. The 
upper limit (3.7) can be rewritten as A|| < AX^ax/'^, 
while the limit (3.10) means L < AYmax/'^ ■ So, for the 
highest energy particles their dispersion along the shock 
can be at most comparable to their mean free path or 
their maximum drift length. Such particles stream along 
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the magnetic field nearly free, advaneing their positions a 
distance ~ rg in each gyroperiod, and they drift across the 
magnetic field changing their y coordinate rg at each 
gyroperiod, as well, provided _B2 3> -Bi- Thus for a given 
minimum acceleration time (3.2) or (3.4) such highest en- 
ergy particles can move comparable distances along and 
across the magnetic field at the shock. 
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